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ABSTRACT 

We have carried out a survey of compact star clusters (apparent size <3") in the southwest part of 
the M31 galaxy, based on the high-resolution Suprime-Cam images (17.5' x 28.5'), covering ~15% of 
the deprojected galaxy disk area. The UBVRI photometry of 285 cluster candidates {V < 20.5 mag) 
was performed using frames of the Local Group Galaxies Survey. The final sample, containing 238 high 
probability star cluster candidates (typical half-light radius rh ~ 1.5 pc), was selected by specifying a 
lower limit of > 0.15" (> 0.6 pc). We derived cluster parameters based on the photometric data 
and multiband images by employing simple stellar population models. The clusters have a wide range 
of ages from ~5 Myr (young objects associated with 24 ^m and/or Ha emission) to ~10 Gyr (globular 
cluster candidates), and possess mass in a range of 3.0 < log(m/m0) < 4.3 peaking at m ~ 4000 m©. 
Typical age of these intermediate-mass clusters is in the range of 30 Myr ^ t < 3 Gyr, with a prominent 
peak at ~70Myr. These findings suggest a rich intermediate-mass star cluster population in M31, 
which appears to be scarce in the Milky Way galaxy. 

Subject headings: galaxies: individual (M31) - galaxies: star clusters - stars: formation 



1. INTRODUCTION 

The Andromeda galaxy, M31, as the nearest spiral 
galaxy, has been regarded to be the most suitable to 
provide supplementary data to the Galactic data for un- 
derstanding the galaxy structures, and to further eluci- 
date the evolution history of disk galaxies. Since a star 
cluster normally represents a single population, many de- 
tailed studies of star clusters in M31 were devoted to de- 
riving information related to the structure and evolution 
of star c l usters and of hosting galaxy: se e the review in 
iKodairal (|2002f ) and the introduction of I Caldwell et al.l 
l|2OO90 . and the works cited therein, among others. 

Recently, many observational works w ere devoted to 
study globu lar clusters (GCs; see, e.g., iBarmby et al.l 
200a [2007L and refere nces therein). The study by 



Cohen et al.l (|2005D . based on high-resolution Keck imag- 



ing, however , suggests that M31 cluster samples (e.g., 
iGalleti et al.l [2007) might be heavily contaminated by 
ast erisms. Never theless, the discussion on this issue 
by I Caldwell et al.l (|2009) states that appearance of low- 
mass clusters dominated by several bright stars depends 
strongly on the photometric passband used, and empha- 
sizes the need of multiband observations to identify clus- 
ters correctly. 

M31 disk cluster studies benefited from the high- 
resolution Hub ble Space Telescope ( HST) images, which 
were used by iWilliams fc Hoded (|2001al lbl) to identify 
young star cluster candidates. Recently an extensive 
archival BVI image survey of HST WFPC2/ACS in- 
strument fields, which are scattering over a wide area of 
M31 in patches, was conducted in search of new clus- 

^ Institute of Physics, Savanoriy 231, Vilnius LT-02300, Lithua- 
nia; donatas.narbutis@fT.vu.lt 

^ Vilnius University Observatory, Ciurlionio 29, Vilnius LT- 
03100, Lithuania 

3 The Graduate University for Advanced Studies (SOKENDAI), 
Shonan Village, Hayama, Kanagawa 240-0193, Japan 

Author to whom any correspondence should be addressed. 



ters (see iKrienke fc Hoded 120071 120081 and references 
therein). 

However, until the recent years, when high-resolution 
wide-field images became available, much less efforts have 
been devoted for a homogeneous wide-field photometric 
survey of objects, which have properties in-between of 
globul ar and open clusters. Our survey, using Suprime- 
Cam (jMivazaki et al.l [20021 ) mosaic images, completely 
covers the area of 17.5' x 2 8.5' in the sout hwestern part 
of the M31 disk (see Kodaira et all 120041 ). A dozen of 
HST frames are located in our field covering less than 
20%. 

Our initial study of M31 compact clusters of an appar- 
ent size less than ^^3" using Suprime-Cam revealed that 
some of the previously suspected GC candidates should 
be classified as open clusters (OCs), suggesting that there 
is a new kind of star cluster, which is not well known in 
the Milky Way galaxy (MW) , but whose counterparts are 
probably present in Large Magellanic Cloud/Small Mag- 
ellanic Cloud (LMC/SMC; Kodaira 2 00l iKodaira et al.l 
I2OOI [2008h . They are apparently more massiv e than 
typic al OCs in the solar neighborhood (Piskunov et al.l 
"200 1, but less massive than the t ypical GCs in MW 
(McLaughlin fc van der Marell [2005[ ) . This suggests the 
presence of an intermediate disk population of star clus- 
ters, a detailed study of which may contribute to the 
understanding of disk evolution history. For examples 
of the typica l images of compact clusters, see the atlas 
pubhshed bv IKodaira et al.l (|2004f l. 

About 50 bright (17.0 S> V 5, 19.0 mag) me mbers 
of our M31 cluster sample (jNarbutis et al.l 120061 ) were 
studied in detail concerning their radial sur face b right- 
ness profiles using a tidal-cutoff King model (King 196^ 
and a power-law EFF model (Els on et a l. 1987 ). Their 
distribution in the half-light r adius versus the absolute 
V^-band magnitude diagram (jSableviciute et al.l 120071) 
shows that these compact cluste rs are occupying th e 
same domain as faint GCs (e.g., iBarmbv et al.l [20071 ). 
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A comparison in the parameter domain of the King 
model (jSableviciute et al] l2006f ) confirms a visual im- 
pression that compact clusters have similar radial pro- 
files to GCs, except for minor irregularities, which be- 
come visible in high- resolution images. Comparison 
of the EFF model param eters of compact clusters and 
massive clusters in LMC (iMackev fc Gilmorell2003D and 
SMC ()Hill fc Zaritskvl l2006f) shows that the present sam- 
ple of compact clusters may have the same structural 
nature as massive clusters in the Magellanic Clouds 
(|Sableviciute et al.|[200l) . 

These promising results motivated us to extend the 
survey of clusters to a deeper magnitude oi V < 
20.5 mag. The homogeneous UBVRI photometric data 
and multiband images for 285 selected cluster candidates 
with absolute magnitudes down to My < —4 m a g wer e 
analyzed in the catalog paper by iNarbutis et al.l (|2008f ) . 
Here, we present the results of simple stellar popula- 
tion (SSP) model fitting using the multiband photom- 
etry data, giving mass and age of these objects in the 
surveyed field of M31. 

We have developed and applied a method to derive 
structural parameters from the surface brightness distri- 
bution of clusters by properly accounting for the contam- 
inating stars, superposed on their smooth profile. The 
derived half-light radii were studied in connection with 
the SSP model fitting results. We have also modeled 
stochastic effects of bright stars (iDeveikis et al.l 120081) 
and briefly discuss their influence on the derived star 
cluster parameters. 

In the following, we describe observations and sam- 
ple selection in Section [2l present structural and evolu- 
tionary star cluster parameters in Section [3l and analyze 
properties of studied sample of M31 cluster candidates 
in Section ID Properties of a few particular objects are 
discussed in the Appendix. 

2. OBSERVATION 
2.1. Data 

We have identified high probability star cluster can- 
didates (hereafter star clusters) by visually inspecting 
the Suprime-Cam y-band mosaic image made from 5 x 
2 minute individual exposures of the 17.5' x 28.5' field, 
centered at ajaooo = 0*^40.9™, Jjaooo = +40°45'. The 
characteristic full width at half maximum (FWHM) of 
stellar images (point-spread function, PSF) is 0.7". See 
iKodaira et al.l (|2004l ) for a description of the Suprime- 
Cam observations and data reduction details. 

We have produced the homogenized UBVRI pho- 
tometric catalog of 285 selected star clusters {V < 
20.5 mag) used as an initial sample in this study, based 
on photometry performed on t he Local Group Ga laxies 
Survey (LGGS) mosaic images (|Massev et al.l[2006l ). The 
variable PSFs of four overlapping LGGS field (F6-F9) 
mosaic images were homogenized to the same FWHM = 
1.5". The mosaics were photometrically calibrated by us- 
ing individual CCD color equation s and the stellar pho- 
tometry catalog from lMassev et all (pOOd) as a local pho- 
tometric standard for zero point calibration. Different 
aperture sizes and photometric background estimation 
areas were selected for ind ividual clusters, as t hey are 
given in the catalog paper (|Narbutis et al.ll2008f ). 

We note that a comparison of three published stellar 



photo metry data sets in our survey field (jNarbutis et al.l 
l2006f ) makes one to be cautious when using tertiary 
standards as local photometric standards. However, 
we used a carefully reduced and calibrated data set by 
IMassev et al.l ()2006D . which passed an internal consis- 
tency check, making it the most accurately calibrated 
photometry catalog of the M31 galaxy to date. 

The color images of clusters were produced from the 
LGGS (U,B,V,I,Ha), GALEX (FVY, NUV), 2MASS 
(J, H, Kg), Spitzer (24 /im), and HI (21 cm) images. They 
served as references for object properties and environ- 
ment study. Seventy-seven star clusters from our catalog 
are located in the HST archive frames. For object se- 
lection, photometry, and multiband image details, see 
INarbutis et"al] (|2008f) . 

We used Suprime-Cam F-band mosaic image 
(FWHMpsF ~ 0.7", image scale 0.2"/pixel) to derive 
structural parameters of star clusters. The LGGS 
field's F7 F-band mosaic image (FWHMpsf ~ 0.8", 
image scale of 0.27"/pixel) was substituted for objects 
saturated or having defects in the Suprime-Cam frame. 
Since the variability of FWHMpgp was only ~5% across 
the mosaic (see Section 13.11 for details), we have applied 
the DAOPHOT package (St etswUfTOSTf) from the IRAF 
program system (|Todv..l993i) to compute a single mosaic 
PSF. IRAF's program "seepsf" was used to create PSFs 
suitable for our structural parameter fitting program. 

The following M31 parameters are adopted in this 
study: distance modulus of m — M = 24 .47 (785 kpc, 
thus 1" = 3.8pc: .McConnachie et al.ll2005f ): center coor- 
dinates aj2ooo = 0^42™44.3% (5j2ooo = 41°16'09" (NASA 
Extragalactic Database); major axis position angle of 38° 
0c Vaucoulcurs 1958): disk inclination angle to the line 
of sight of 75° (iGordon et al.ll2006D . 

2.2. Sample Selection 

The star clusters were selected from the initial sample 
of 285 objects by considering the following criteria: (1) 
the derived half-light radius is of rji > 0.15" (> 0.6 pc) — 
this limit arises due to the resolution of ground-based ob- 
servations; (2) there are no nearby contaminants strongly 
influencing the accuracy of the derived structural and 
evolutionary parameters. Note however, that several 
genuine star clusters, judged from the HST images, did 
not meet these criteria. 

A few dim, i.e., low central surface brightness, ob- 
jects with an unreasonably large derived half-light ra- 
dius or composed of several resolved stars were rejected 
because of low reliability of determined parameters. Al- 
though PEGASE (|Fioc fc Rocca-Volmeranga[l99l SSP 
models incorporate nebular emission lines for young age 
objects, we stress that structural and evolutionary pa- 
rameters of clusters showing up Ha emissions or super- 
posed on inhomogeneous diffuse Ha background should 
be considered with caution. Note however, that objects 
with strong Ha emi ssion were cataloged separately by 
IKodaira et al.l (j2004l ). therefore, these objects were ex- 
cluded from the analysis in this study. We note, that 
foreground starburst galaxies possessing 24 fim emission, 
which have been reve aled during spectroscopic study by 
I Caldwell et al.l ()2009f ). might be present in a sample of 
young {t < 20 Myr) clusters (see Appendix for details on 
KW271 and KW279). 
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This cleaning yielded a sample of 238 star clus- 
ters. Their ID num bers and center coordinates from 
iNarbutis et alj (120081 ) are supplemented with the struc- 
tural and evolutionary parameters described in the fol- 
lowing, and are presented in Table 1, the full form of 
which is available in the electronic edition of the Jour- 
nal. 

3. MODEL FITTING 
3.1. Structural Parameters 

The studied star clusters, projected on the crowded 
disk of M31, are unresolved or semiresolved in 
the Suprime-Cam images. Background and/or clus- 
ters' bright stars superposed on clusters' surface 
brightness distribution strongly influence photometric 
(jNarbutis et al.l r2007bl ) and structural parameter accu- 
racy. However, stars projecting nearby to the cluster 
cannot be simply masked out because of flux superpo- 
sition in their extended wings. To avoid this, we have 
developed a method to derive structural parameters of 
star clusters by including resolved stars in the cluster 
model fitting procedure. The method is presented in D. 
Semionov et al. (2009, in preparation); here we provide 
its brief description and derived half-light radii of star 
clusters (Table 1). 

The developed program tool models star clusters as a 
smooth analytical surface brightness distributions. Ob- 
servational effects are taken into account by convolving 
cluster models with PSF. Individual resolved stars are 
modeled with PSF. Free parameters, which describe a 
shape and position of the cluster model, are fitted using 
a genetic algorithm, providing the global solution. Fine 
tuning of this solution is performed using Levenberg- 
Marquardt nonlinear least-square algorithm. A photo- 
metric background value, which is kept constant during 
model fit, is estimated within individually selected cir- 
cular annulus centered on the object. The model fitting 
radius is chosen seeking to enclose extended wings of the 
surface brightness distribution. The residual images are 
analyzed to ensure that there are no visible systematic 
deviations. 

We have fitt ed 285 clusters from the initial sample with 
circular King pa^[l96l ) and EFF (|Elson et all 119871) 
models. Since these models have different sensitivity to 
the cluster's "core" and "wings," comparison of the half- 
light radii derived by their means provides the robust- 
ness of the model fit and r/j. Elliptical models were not 
used in this study due to stochastic distribution of bright 
cluster stars, whic h introduce fake ellipt icity for the low- 
mass clusters (see iDeveikis et al.l 120081 for a discussion 
of stochastic effects). A number of contaminating stars 
used for the model fitting was: in 40%, 1 in 30%, and 2 
in 20% cases. The maximum number of fitted stars was 
6. 

Typical fitted values of structural parameters are: (1) 
the core radius, rc 0.75 pc, the concentration param- 
eter, c ~ 1.5 (the King model); (2) the effective radius, 
~ 0.75 pc, and the power-law index, n ~ 1.3 (the EFF 
model). For ~70% of clusters the derived c and n pa- 
rameter values were confined to the reasonable li mits of 
0.5 < , c < 2.5 (see boundaries for MW GCs in iHarrisI 
Il996f ) and 1.05 < n < 3.5, respectively. For practical 
convenience, the remaining ~30% of clusters, having de- 



termined parameters outside these limits, were re-fitted 
by setting particular parameters to the corresponding 
boundary values. We note however, that although this 
slightly affects derived Tc or rg, it does not change the 
derived half-light radii nor introduces systematic devia- 
tions in the residual images. 

Finally, intrinsic half-light radii of objects, based on 
the King, r^, and the EFF, rf, models were derived. 
They are plotted in Figure 1(a) and show a good corre- 
lation. The average half-light radii values, rh = {r^ + 
rf )/2, for 238 star clusters are provided in Table 1. 

To estimate both the robustness of the model fitting 
and the uncertainty of the derived half-light radius, cr^^ , 
we have performed two tests: (1) varied photometric 
background value, /isky, which is kept constant during 
model fitting; and (2) decreased fitting radius, rat, from 
the adopted initial value. These tests indicate that r^ 
based on the King and the EFF models are in a good 
agreement, and the overall accuracy of the rh value is 
better than 20%. The typical is indicated by error- 
bars in Figure 1(a). Based on these results, we ap- 
plied a conservative lower limit (which is mainly a sub- 
ject to the resolution of ground-based observations) of 
rii — 0.15" (0.6 pc) for the sample selection. Note how- 
ever, that some clusters below this limit have been iden- 
tified in the HST images, which turned out to be of a 
very compact nature or dominated by a central bright 
star. 

To test the influence of PSF variability on the de- 
rived half-light radii, we have divided survey fleld into 
15 regions and constructed individual PSFs for these re- 
gions. These PSF models were treated as "virtual ob- 
jects" and analyzed as star clusters in attempt to derive 
their structural parameters. The derived half-light radii 
of eight "virtual objects" are displayed with asterisks in 
Figure 1(a) and are confined in the range of ru < 0.4 pc. 
Note that only eight of 15 "virtual objects" are visible, 
since others failed to be fitted due to FWHM smaller 
than that of the adopted mosaic PSF. The scatter of 
rh < 0.4 pc corresponds to the rh uncertainty of 20% for 
objects with rh ^ 2 pc. It is the upper limit of uncer- 
tainty arising due to PSF variability across the survey 
field. For some individual clusters, we have additionally 
derived structural parameters by using local and mosaic 
PSFs, and found results being in a good agreement. 

In case of semiresolved star clusters, subtraction of con- 
taminating stars has a strong effect on the accuracy of 
derived half-light radii. We have automatically included 
all stars within rgt -I- FWHMpsF from the cluster's cen- 
ter and having magnitude < (mdustor + 3) mag, i.e., up to 
3 mag fainter than the cluster itself, into the fitting pro- 
cedure. Also some individual strong contaminants were 
marked by hand. Conventional tests of varying the back- 
ground level and the fitting radius indicate that ^--^20% 
accuracy of rh is achieved. Visual inspection of resid- 
ual images proved that clusters and contaminating stars 
are subtracted well. We have estimated the goodness of 
model fitting by comparing the standard deviation of the 
sky background within the cluster's area after model sub- 
traction and that of the nearby sky region. For ~90% of 
cases, these values were equal, except for the cases with 
larger residuals visible in the cluster's core, which, most 
probably, arise due to semiresolved stars. They are the 



Vansevicius et al. 



[pc] 




Fig. 1. — Half-light radii of star cluster candidates, (a) De- 
termined half-light radii of 285 initial sample objects — based on 
the EFF model, r^, and the King model, rj^ (squares mark re- 
jected cluster candidates; circles — selected 238 clusters; asterisks — 
PSFs constructed in the survey field and analyzed as "virtual ob- 
jects;" dashed lines — the applied half-light radii selection limit 
r/i > 0.6 pc; error-bars — typical cr^^). (b) Histogram of average 
half-light radii, = {r^ +r^)/2, of 238 clusters, overplotted with 
a slope r] = —2.2 line, defined as N{rf^)drf^ oc r^dr^. 

main source of uncerta inty for low-mass yo ung star 
clusters (see Figure 3(c) in lDeveikis et al.|[2008D . 

We have compared rh values of 50 bright objects from 
our sa mple with the ones derived by ISableviciiite et aTl 
(j2007[) . Despite different analysis methods employed, we 
have found a good agreement between these two studies. 
The maximum difference of ~30% is attributed to the 
influence of contaminating stars not properly accounted 
for in the previous study. 

3.2. Evolutionary Parameters 

The UBVRI aperture CCD photometry data of 285 
star clusters from the initial sample (Narbuti s et al.l 

20081) were compared with SSP models (PEGASE code: 
Fioc fc Rocca-Volmerangelll997D assuming the standard 
initial mass function (IMF: iKroupal [2003 ) and various 
metallicities to yield the best fit model evolutionary pa- 
rameters, i.e., simple SSP model fitting was performed. 
The interstellar extinction w as accounted for ass uming 
the standard extinction law (jCardelH et al.|[l989l ). The 
adopted parameter quantification technique and intrin- 
sic precision analysis, based on the UB VRI integrated 
photometry, was presented in detail by iNarbutis et al] 
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Fig. 2. — Accuracy of the derived evolutionary parameters of 
238 star clusters. Panels (a) and (b) show the standard deviation 
of age, CTiog^i/jyiyr) 1 plottcd VS. V-band magnitude and vs. age, 
log(t/Myr), respectively. Panels (c) and (d) display the histograms 
of standard deviations: color excess, i^e{B—V)^ ^.nd metallicity, 
""[M/H]) respectively. 



(|2007aD and lBridzius et all (|2008D . This technique, sup- 
plemented with a detailed investigation of multiband im- 
ages, introduced to reduce the age-metallicity-extinction 
degeneracy, was used in the study. 

The simple SSP model fitting results for 238 star clus- 
ters are provided in Table 1. The determined evo- 
lutionary parameters are: absolute T^-band magnitude 
(corrected for aperture correction, see below for de- 
tails), My; age in Myr, log(f/Myr); mass in solar-mass 
units, log(TO/m0); metallicity, [M/H]; and color excess, 
E{B - V). 

The standard deviation of age, cf\a^i^t/Myc)^ is plotted 
versus F-band magnitude and age, log(t/Myr), in Fig- 
ures 2(a) & (b), respectively. The standard deviation 
of age increases typically from '^0.1 dex for young to 
~0.3dex for old objects. Wc note that systematic dif- 
ferences of the derived parameters could be expected if 
another SSP model bank would be used. The histograms 
of the standard deviation of color excess, (Je{b-v)^ ^-nd 
metallicity, (T[m/h], are shown in Figures 2(c) & (d), re- 
spectively. The characteristic cte{b-v) is ^0.05, reach- 
ing ~0.15 for t ~ lOGyr objects due to stronger age- 
extinction degeneracy. The typical (Tjm/h] is ~0.3dex, 
reaching up to '--'0.7 dex for t ~ 200 Myr objects due to 
strong er age-metallicity degeneracy (see .Narbutis et al.l 
l2007al for parameter degeneracy maps). 

We note that in this study numerous SSP model sim- 
plifications were assumed: (1) a constant IMF; (2) a 
standard for MW and constant for all cluster stars ex- 
tinction law; (3) solar element ratio. Additionally, in 
cases of unreliable photometric solutions for metallicity, 
we restricted [M/H] determination within lower and up- 
per boundaries as a function of age taken from the chem- 
ical evolution model p rediction . Referring to models of 
MW and M31 bv iRenda et all (|2005i") . we assumed that 
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log(t/Myr) log(m/mQ) 

Fig. 3. — Parameters of 238 star clusters, (a) Half-light radius, 
Vf^, plotted vs. age, log(i/Myr). (b) rj, vs. mass, log(m/mQ), 
overplotted with lines of a constant half-mass density, ph = 1, 10, 
10^, 10"^, and 10^ tuq pc~^, from left to right, respectively, (c) and 
(d) Histograms of the derived cluster age, log(t/Myr), and mass, 
log{m/m0), respectively. In panel (b) error- bars indicate typical 

metallicity evolution in the M31 survey field is not signif- 
icantly different fro m that of solar neighborhoo d, which 
has been modeled bv lSchonrich &: BinnevI (|2009( ). There- 
fore, the [M/H] boundaries roughly between —1 and -1-0.5 
were considered, also taking into account possibility that 
oldest objects could be globular clusters in M31 halo with 
[M/H] as low as -2. 

The mass of clusters in the solar-mass unit, m/iriQ, 
was calculated equating the mass-to-luminosity ratio of 
the best fit SSP model of age, t, and metallicity, [M/H], 
to the absolute V-band magnitude, My. The mass of 
clusters strongly depends on the derived age, therefore, 
the primary sources of mass uncertainty are the accuracy 
of age and extinction. Since cluster candidates in the 
crowded field were measured through small apertures, we 
have taken into account an individual aperture correction 
for each star cluster by applying the following procedure. 
Artificial clusters based on the best-fitted King model pa- 
rameters were generated and convolved with the homog- 
enized PSF of FWHM = 1.5"; i.e., resolution of model 
cluster images was matched to the LGGS frame resolu- 
tion employed for cluster photometry. Individual aper- 
ture used for a real cluster was applied to measure corre- 
sponding model cluster, the difference between measured 
and total flux was assumed as an aperture correction for 
a real cluster photometry. A typical l^-band aperture 
correction is of '^0.4 mag and translates into a correction 
of mass by '^0.15 dex in log(m/mQ) scale. 

Analysis of star cluster population requires a proper as- 
sessment of selection effects and res ulting completen ess. 
Stochastic modeling of clusters ( De veikis et al.ll2008f ) un- 
folded a broad variety of their possible appearances, indi- 
cating that visual selection of low-mass clusters could be 
biased. However, for older clusters, i.e., after the phase 
of supergiants {t > 100 Myr), this problem appears to 



become less crucial. 

4. RESULTS AND DISCUSSION 

Here, we analyze the structural and evolutionary pa- 
rameters of 238 star clusters in the southwestern field of 
the M31: half-light radius; age; mass; color excess; metal- 
licity; and spatial distribution. Prior to presenting our 
results, we stress that important sources of systematic 
bias and uncertainty of the derived parameters, in par- 
ticular for young low-mass star clusters, are stochastic 
effects (Deveikis et al. 2008). Assuming stochastic na- 
ture of star formation to randomly populate the stellar 
IMF, we modeled influence of stochastic effects on mea- 
surable cluster characteristics, and found that structural 
and evolutionary parameters could be biased — for some 
cluster age and mass i nterv als — in a systematic way. Re- 
cently, [Barker et al.l ()2008r ) gave a strong caution that 
the standard SSP model analysis significantly underesti- 
mates uncertainty in the derived cluster age. However, 
only standard SSP model fitting was performed in this 
study; see a recent attempt to solve this proble i n for star 
clusters of solar metallicity bv lMaiz Apellani"3 (|2009f ). 

4.1. Half-Light Radius 

The histogram of clusters' half-light radius, r^, dis- 
played in Figure 1(b), spans the range from ~0.6pc to 
^10 pc with one large object of rh ^ 14 pc (KW249), 
which is described in the Appendix; the distribution 
peaks at ^ 1.5 pc. For comparison, the peak of the 
half-li ght radius distribution of MW GCs from iHarrid 
(fToM February 2003 rev.°) catalog is at ~2. 5 pc. There- 
fore, the majority of clusters from our M31 sample are 
smaller (more compact) than typical MW GCs. 

The turnover and a lack of small-size clusters close to 
the applied lower limit of r/i = 0.6 pc can be attributed to 
the visual selection of the cluster sample. The exponen- 
tial decrease in the number of clusters in the larger half- 
light radii domain should be real. Note however , tests 
performed with "SimClust" (jPeveikis et al.l[2008| ) show 
that only extended clusters of low luminosity (My ^ 
—4.5) projected on star-forming regions could be missed 
by visual inspection. The overplotted power-law slope 
77 = —2.2 line, defined as N{rfi)drh ex. r^d rh, was found 
as the best fit for M51 star clusters by iBastian et al.l 
HqM), and agrees well with our star cluster size distri- 
bution in the range oirh = 3 — 8pc. However, a consider- 
ab ly steeper slope was found in the M51 cluster sample 
of iScheepmaker et all (l2007f). which is five times more 
numerous than that of IBastian et al.l (|2005[ ). 

The half-light radius of clusters is plotted versus age, 
log(t/Myr), in Figure 3(a). No obvious half-light radius 
evolution with the age o f cluster population, as has been 
found in LMC/SMC bv iMackev fc Gilmord pOOl . can 
be judged from this diagram, since objects with large rh 
are observed over a wide range of age. However, the lack 
of small-size clusters at age t ~ 3 G yr can be a ttributed 
to the selection effects discussed bv lHunter et al. (200^ 
for the LMC/SMC case, such that with increasing age 
only bright and, therefore, massive extended clusters are 
selected in the magnitude limited sample. At i ~ 10 Gyr, 
halo GCs likely dominate in our sample and show 
distribution similar to that of the MW GCs population 

^ See jhttp:// phy swww. mcmaster . ca/ ~harris /mwgc . dat [ 
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with rfi as small as ^^Ipc. They have high densities (see 
Figure 3(b)) and likely have survived tidal disturbances 
by the disk. 

The half-light radius versus mass, log(TO/m0), distri- 
bution is shown in Figure 3(b), overplotted with five lines 
of constant half-mass density, ph = 3m/87rr)^, ranging 
from 1 to IO^ttiqPC"^. The most numerous clusters in 
the M31 southwestern field are compact clusters of small 
size r/i ~ 1.5 pc, intermediate mass m ^ 4000 m©, and of 
young t ^ 100 Myr age. Largest clusters seem to be en- 
veloped by the fines of ph from 1 to 100 pc~^, lending 
an impression that massive clusters tend to be bigger in 
size as far as clusters of relatively low specific density of 
Ph < lOOm0pc~'^ are concerned ([Hunter et al.ll2003l i. 

Comparing the versus mass distribution of 

our sample with a sa mple of M51 star clusters 
(jScheepmaker et aLll2007D . we note a lack of objects with 
4 < log(m/m0) < 5 and r/j < 5pc in our sample. The 
most massive, log(m/m0) > 4.5, clusters are "branch- 
ing" in Figure 3(b) at small size, < 3pc, and high 
density, ph ^ 10^ pc~^, which is characteristic to old 
GCs. 

The oldest clusters, which have highest densities, are 
most likely halo GCs of small r^,. That the typical life- 
time of 10** TTiQ mass star cluster is estimated to be of 
^300 Myr in the survey field (see Section [4771 for details) 
implies very low probability for disk star clusters to reach 
an age of ~10Gyr without disruption. However, radial 
velocity measurements are necessary to confirm whether 
these objects belong to halo. 

4.2. Age 

The distribution of clusters in the mass versus age di- 
agram is shown in Figure 4. The upper right domain of 
the diagram, t > 3Gyr and 4.0 < log(m/m0) < 7.0, is 
occupied by ~30 classical GC candidates. The remain- 
ing ~210 are most likely disk clusters, which reside in the 
domain oi t < 3Gyr and log(TO/mQ) < 4.5. The curve 
along the lower envelope of cluster distribution is due to 
the selection limit of < 20.5 mag, which was calculated 
with PEGASE SSP models of [M/H] = -0.4 dex metal- 
licity, assuming zero extinction. Note, that star clusters 
are located above the selection limit due to aperture cor- 
rections taken into account for mass calculation. How- 
ever, the upper envelope of the distribution should be 
fr ee of selection effects. 

iCervino fc Luridianal ([2004f ) have shown that SSP 
models cannot be used in a straightforward determinis- 
tic way for small-mass cluster analysis, when integrated 
luminosity of a model is lower than that of the brightest 
star included in the used isochrone, i.e., "Lowest Lu- 
minosity Limit" (LLL) requirement. They defined the 
smallest mass, m}-'^^, associated with the LLL, and the 
mass, m'^^^, corresponding to the 10% stochastic fluctu- 
ation of the 1^-band luminosity, c^L^, = 0.1 x Ly, which 
translates into ~0.04dex uncertainty of log(m/m0). We 
have overplotted both_parameters as a function of age in 
Figure 4, taken from iCervino fc Lurid iana (2004,), noting 
that m'^^^(t) ~ 10 X m^^^{t). No te however, that SSP 
models of solar metallicity used bv ICervino fc Luridianal 
((2004) do not significantly influence a qualitative assess- 
ment of mass versus age distribution. 

Judging from Figure 4, m^^^ is log(m/m0) ~ 4 for 
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Fig. 4. — Mass vs. age diagram of 238 star clusters, overplotted 
with sample selection limit, V < 20.5 mag, (solid line), the Low- 
est Luminosity Limit, m^^^(t), (dashed line), and 10% stochastic 
fluctuation of cluster V-band luminosity, m°'^^{t), (dotted line). 

clusters younger than '^lOMyr. At the age of '^50Myr, 
^j^LLL (jj-Qpg log(m/TO0) ^ 3, therefore, the majority 
of clusters from our sample satisfy the LLL requirement. 
At t 100 Myr, the m^^^ coincides with V < 20.5 mag 
selection limit; older objects are well above LLL. There- 
fore, clusters with log(m/m0) > 4.5 should have the 
most accurate evolutionary parameters. 

The stochastic effects on the derived cluster mass can 
be estimated by referring to the dotted line in Figure 4, 
based on the Ccrvino & Luridiana (2004) methods. The 
line displays 10% uncertainty of cluster's luminosity at a 
specific mass for a given age. Since we derived the mass 
from the y-band luminosity, the mass as a function of 
age, having 10% uncertainty of mass, is depicted by the 
dotted line in Figure 4. Formally, mass uncertainty in- 
creases to 100% at the LLL indicated by the dashed line. 
However, the uncertainty of mass is additionally affected 
by the uncertainty of the derived age, which can be eval- 
uated from the slope of solid line in Figure 4, and also by 
the uncertainty of extinction shown in Figure 2. There- 
fore, the LLL provides only the lower limit for uncer- 
tainty arising due to stochastic effects, which primarily 
bias age and extinction, and consequently — mass, i.e., 
actual accuracy of evolutionary parameters is lower than 
it is shown in Figure 2. 

The histogram of cluster age, displayed in Figure 3(c), 
spans the range from t ^ 5 Myr to ~12 Gyr. Since evolu- 
tionary parameters of the clusters with Ha emission are 
less reliable due to emission lines altering photometry 
and stochastic effects, we cannot properly estimate pa- 
rameters of the clusters younger than ~20 Myr. There- 
fore, we did not specifically target our survey for such 
star clusters. 

The age distribution of clusters is displayed in Fig- 
ure 5(a). We have calculated the number of clus- 
ters per linear age interval, dN/dt, by binning data in 
dlog(i/Myr) = 0.2, 0.3, and 0.5 age intervals, while also 
shifting them by a half-bin width. These results were 
averaged and smoothed. The symbol size indicates an 
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Fig. 5. — Age and mass distributions of 238 star clusters, in- 
dicated by open circles (size indicates tlie uncertainty of data 
binning) in both panels, (a) Number of clusters per age inter- 
val, log(dA''/dt), vs. age, log(t/Myr) (circles), overplotted with 
power-law slopes expected for cluster population forming continu- 
ously and: (1) fading evolutionary below the detection limit (solid 
line); (2) instantaneously disrupting at some arbitrary time (dash- 
dotted line), (b) Number of clusters per mass interval, normalized 
to the total number of clusters in the sample, log(dA'^/dr)i/A'^tot ), 
vs. mass, log(m/m (7)) (circles), overplotted with Schechter's func- 
tion (dashed line; IGielesI 120 09). adopting the turnover mass of 
m* = 2 X 10^ m0. Stars correspond to the intermediate age star 
clusters from our sample in-between of ^^100 Myr and ~3 Gyr; filled 
circles cor respond to the star c lusters from MMT spectroscopic 
survey by ICaldwell etU] l|2009l . subsample V < 19.0 mag, ages 
in-between of ~ 100 Myr and an upper age limit of their sample 
~lGyr). 

uncertainty of the resulting distribution, which is plot- 
ted versus age in a log scale. The age distribution is 
reasonably well represented with two power-law lines of 
different slopes, and is to be interpreted as a result of 
star cluster formation/disruption processes (see Section 
14.71 for discussion). 

The cluster number distribution dN/dt (Figure 5(a)) 
in general follows the "evolutionary fading" line. How- 
ever, after subtracting the "evolutionary fading" number 
gradient, a gradual increase in the number of clusters per 
log age interval from t ~ 20 Myr to t ^ 100 Myr is ob- 
vious. A peak of the cluster number distribution is fol- 
lowed with a decrease in the number of clusters, which 
is a consequence of an evolutionary cluster fading effect, 
being always significant in the magnitude limited sample 
(see Figure 4), and star cluster disruption. We note that 
analysis of both Figures 3(c) & 5(a) points to the peak 
of the cluster age distribution at i ~ 70 Myr, suggesting 
an enhanced cluster formation episode at that epoch. A 
smaller secondary peak at ^10 Gyr is due to old GC can- 
didates. 

4.3. Mass 

The histogram of cluster mass, displayed in Fig- 
ure 3(d), spans the range from ~10^ to ~5 x 10^ m© and 
shows a prominent peak at m ~ 4000 TOq. The decrease 
in the number of clusters in the low-mass domain is ap- 
parently not physical, but arises due to selection effects. 
The estimated complete ness of our cluster sa mple is of 
~50% at F ~ 20.5 mag ()Narbutis et al.ll200l . Judging 
from these circumstances, the completeness at the age of 
~100Myr and the mass of ~4000rno (Figure 4) should 
be credible, even taking into account a stochastic scatter- 
ing of cluster luminosity at this age and mass. Therefore, 
decline in the number of clusters per log mass interval in 
the high-mass domain should be real. 



The mass distribution of clusters is displayed in Fig- 
ure 5(b). The number of clusters per linear mass interval, 
normalized to the total number of A^tot = 238 clusters in 
our sample, dN/dm/Ntot, was calculated by binning data 
in dlog(m/TO0) — 0.2, 0.3, and 0.5 mass intervals, and 
also shifting them by a half-bin width. These results were 
averaged and smoothed. The symbol size indicates the 
uncertainty of the resulting distribution. For a more de- 
tailed comparison of intermediate age (from ^100 Myr to 
~3 Gyr) star clusters from our sample we have used the 
selected subsample {V < 19.0 mag, age from ~100Myr to 
an upper age limit of their sa mple ^1 Gyr) of star cl usters 
from the M31 disk study bv ICaldwell et all (I2009D The 
dN/dm distribution presented in ICaldwell et all (pool ) 
subsample was multiplied by 0.15, taking into account 
that our survey field covers only ~15% of the M31 disk. 
We note, that there are only e ight objects in common 
between our sample and that of ICaldwell et al.l (120091. 

De cline in the distribution of mass (jCaldwell et al.l 
|2009() less than log(m/TO0) ~ 4.3 (see Figure 5(b)) arises 
due to a sample incompleteness and resembles a shape 
of our intermediate age sample distribution for mass less 
than log(m/TO0) 3.5. Although we have only few star 
clusters in the high-mass range, th e good match with a 
dN/dm slope (jCaldwell et al.ll2()09l ) and extension down 
to the mass as low as log(m/m0) ~ 3.7 allows us to con- 
strain a shape of the star cluster mass function for age 
in-between of '^100 Myr and ~ 1 Gyr. We o verplotted 
Schechter's mass function from iGielesI (|2009l Equation 
(1)), adopting the turnover mass of m* = 2 x 10^ m©, 
which shows a good a greement with both o ur inter- 
mediate age sample and ICaldwell et aP ()2009| ) subsam- 
ple data. We note that only vertical adjustment for 
the Schechter's function was applied to match data at 
m ~ 10^ niQ. 

Although we lack low-mass star clusters in our sample 
due to the selection limit at V ^ 20.5 mag {My ^ —4.0), 
the detailed HST study of s tar clusters in M31 by 
IKrienke fc Hodei (pOOTl I200I has revealed a numer- 
ous population of star cluster candidates as faint as 
V ^ 23 mag (My ^ —1.5). They are spatially well corre- 
lated with brighter counterparts oit 100 Myr from our 
sample, therefore, it is reasonable to assume that they 
could be of a similar age. Thus, their mass should be 
significantly smaller than ~4OOOm0, corresponding to 
the mass rang e of typical MW OCs in the solar neigh- 
borhood ( Piskunov et al.l 120081 ). Although the parame- 
ters of MW star clusters are derived basing on individual 
stars, note that masses of MW clusters are a subject to 
the accuracy of their derived tidal radii, as discussed by 
iPiskunov et aH (pOOl . 

The high-mass domain of clusters from our sample is 
found to be overlapping well with the sample of ^^140 
young clusters, sharing the M31 disk ro tation, studied 
using MMT spectra and HST images bv ICaldwell et al.l 
(|2009f l. who found a mass peak at log(m/TO0) ^ 4, sig- 
nificantly higher than for our cluster sample, probably, 
due to brighter cluster selection limit applied for spec- 
troscopy. The mass range of clusters from our sample 
younger than ^3 Gyr falls in-betwe en of the MW GCs 
(McLaughlin fc van der Marelll2005h and the MW OCs 
(Piskunov et all 120081 ). indicating that we are mainly 
dealing with the disk clusters from the intermediate mass 
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range. The mass distribution of clusters older than 
~3 G yr overlaps well with the mass distr ibution of MW 
GCs (IMcLaughhn fc van der Marel[2005f ). 

Comparing the mass versus age diagram of our 
cluster sample (Fi gure 4) wi th those of LMC/ SMC 
(|Hunter et all 120031 ) and M51 ijBastian et al.ll2005f l. we 
see a prominent lack of massive, log(TO/m0) ~ 5, clusters 
younger than '^SGyr in M31. Although our survey cov- 
ers only ~15% of the deprojected M31 dis k it contains a 
representative part of its star-forming ring ([Gordon et al.l 
HOOe) in the vicinity of NGC 206. We also note that other 
studies of the M31 disk clusters, which covered the whole 
galaxy, did not detect a significant number of clusters 
more massive than \og,(m / m p,) ^ 4.5 in this age range 
either (| Caldwell et al.ll2009l . and references therein). 

A few massive, log(m/m0) 5, objec ts with an esti- 
mated age of t < 1 Gyr were reported by 'Caldwel l et al.l 
(j2009f ). who found that most of their young clusters have 
a low concentration, typical to the MW GCs. Their spa- 
tial distribution is subject to patchy field selection of 
HST high-resolution imaging over the M31 disk. Al- 
though there are only eight objects in common in their 
and our samples, judging from the overlapping mass 
range and the reasonable agreement of the derived age, 
we may tentatively assume that both samples belong to 
the same category of star clusters. 

The general lack of star clusters of log(m/m0) ~ 5 
mass in the M31 disk co uld be due to a low- average star 
formation rate in M31 (|Barmbv et al.l [2bG6'). However, 
an unfavorable inclination angle of M31 could hide star 
clusters emb edded in the M31 disk, mimicking the situa- 
tion in MW (jClark et al.ll2005D — there could be massive 
star clusters hidden by dust clouds and yet to be detected 
in M31. This hypothesis might be supported by the 
molecular cloud mass d istribution in M31 (Niete n et al.l 
I2006t [Rosolowskvll2007D . where star clusters with a mass 
of up to log(m/TO0) ~ 5 are expected to form, if the 
typical star formation efficiency of ^30% in the cores of 
molecular clouds is assumed. 

4.4. Extinction 

A typical color excess of star clusters in our sample 
is E{B -V) 0.25. Clusters younger than --500Myr 
have a range of E{B — V) < 1.0, which is consistent 
with a large extinction of objects residing in the galaxy 
disk, e.g., four high extinction, E{B — V) ~ 1.1, ob- 
jects of t ^ 100 Myr are projecting close to the 24 fim 
emission regions and prominent dust lanes. The range 
of E{B — V) < 0.4 derived for clusters of intermediate 
age, 500 Myr < i < 3 Gyr, might be related to the evolu- 
tionary fading of clusters and the observed narrow mass 
interval in this age range (see Figure 4), therefore, ob- 
jects possessing color excess higher than ~0.4 are miss- 
ing from our sample. We note that E{B — V) versus 
projected distance from the M31 center does not show 
any tight correlation. Moreover, E{B — V) of an individ- 
ual cluster seems to be more strongly dependent on the 
local interstellar matter environment than on the radial 
distance from the M31 center. 

The GC candidates of t > 3 Gyr have color excess over 
the range of E{B — V)< 0.6, consistent with their pres- 
ence in the M31 halo on both sides of the disk. There 
are numerous objects with E{B — V) ^ 0.04, which is 
in a good agreement with MW color excess in the direc- 
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Fig. 6. — Metallicity, [M/H], of 238 star clusters plotted vs. age, 
log(t/Myr), and projected distance from the M31 center, Rmsi, in 
panels (a) and (b), respectively. In panel (a) error-bar indicates 
typical uncertainty of metallicity taken from Figure 2(d). In panel 
(b) clusters are divided into two age groups of t < 3 Gyr (small 
circles) and t > 3 Gyr (large circles), and overplotted with the av- 
erage metallicity trends for disk red giant stars from Worthc v et al.l 
120051 . dashed line) and iBellazzini et al.i (j2003j, dotted line). 

tion of the surve y field in M31, estim ated from the MW 
extinction maps (jSchlegel et al.l ll998V 

Note howeve r , that for eight objects in common with 
iCaldwell et"all (|2009f ) sample, which have E{B - - 
0.25 derived in their study, we determine color excess 
from 0.05 to 0.65. Although Caldwell et al. (2009) de- 
termined extinction by examining the shape of spectra 
continuum, for those cases when late-type stars domi- 
nate they used a mean color excess for young clusters of 
E{B — V) ~ 0.28, which is in agreement with the typ- 
ical value of E{B — V) ^ 0.25 derived in our study for 
presumably the same population of star clusters. 

4.5. Metallicity 

The metallicity, [M/H], of star clusters, plotted versus 
age and projected distance from the M31 center, i?M3i, in 
Figures 6(a) & (b), respectively, spans a range from —2.0 
to -1-0.5 dex. We divided clusters into two age groups of 
t <3 Gyr and t >3 Gyr, indicated by a different symbol 
size in Figure 6(b). The apparent narrowing of the metal- 
licity scatter at an age of ~200 Myr is attributed to the 
SSP model fitting artifact due to a strong age-m ctallicity 
degeneracy in this age domain (jNarbutis et al., 2007ai). 
Although the accuracy of photometrically derived metal- 
licity is low (see Figure 2(d)), the evolutionary trend with 
time is apparent. 

A rather constant average metallicity of [M/H] ^ 
—0.4 dex is observed for clusters younger than ~lGyr, 
however, there is a slight tendency of the average [M/H] 
to decrease with increasing galactocentric distance. The 
radial profile of cluster met allicity matches well those 
of dis k red giant star s from iWorthev et al.l (|2005D and 
Bella zzini et al.l ()2003[ ) shown by lines in Figure 6(b). 

There are 20 clusters in our sample, which have [Fe/H] 
est imates from t h e Wil liam Herschel Telescope spectra 
by iPerrett et al.l ()2002[ ) . By accounting for a-element 
overabundance typical to MW GCs, we see a reason- 
able agreement between photometric and spectroscopic 
metallicity estimates of old, t ~ 10 Gyr, clusters. Two of 
those metal-rich massive clusters (KW221 and KW225), 
presumably similar to GCs, have metallicity of [M/H] ^ 
—0.8 and 0.0 dex, and are found in our survey field 
area closest to the M31 bulge, located at the projected 
distance of i?M3i ~ 5kpc from the galaxy's center. 
Note however, that for several objects of ~100Myr age. 
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Fig. 7.— Distribution of 238 star clusters in the M31 SW field, 
overlaid on Spitzer (24 /im) image. Clusters are divided into three 
age groups: t < 150 Myr (stars), 150 Myr ^ 3Gyr (circles), and 
t > 3 Gyr (asterisks). Elliptical ring segments, indicating projected 
distance from the M31 center, iJM3l =6 — 18 kpc, are marked with 
dashed lines. North is up, east is left. 



the metaUicity estimated by iPerrett et alJ (|2002f l is of 
[Fe/H] ~ — 2.2dex, i.e., significantly lower than a photo- 
metrically estimated [M/H] — 0.4dex. 

The elabo rated M31 evolution models, such as "M31 
model-b" bv lRenda et al . (2005), may well reproduce the 
metallicity gradient and the metaUicity change in time, 
as they are presented in Figure 6. Detailed studies of 
the age-dependent profiles of the disk cluster population 
would provide further constraints on the M31 disk evo- 
lution models. 

4.6. Spatial Distribution 

Positions of star clusters are indicated on the Spitzer 
(24 /im) image of the M31 southwestern field in Fig- 
ure 7. Clusters are divided into three age groups: young 
(t < 150 Myr); intermediate (150 Myr <t< 3 Gyr); and 
old {t > 3 Gyr), shown by different symbols. The dis- 
tribution of young clusters resembles that of warm dust, 
but is shifted outward to larger i?M3i by ^1.5 kpc along 
the galaxy's major axis. The intermediate age clusters 
are smoothly distributed over the whole field. Old age 
clusters (GC candidates) scatter over the area, but show 
a higher concentration toward the M31 bulge. 

Star clusters with ultraviolet {GALEX) or 24 /im 
{Spitzer) emission are mainly distributed along the sites 



of active star formation. At the bottom of the survey 
field, young and intermediate age clusters appear to be 
well mixed. Around the position of [6', 16'], where a lack 
of warm dust is attributed to the "split" of the star form- 
ing ring, intermediate age clusters dominate. However, 
around [17', 11'] mainly j/oMn^ clusters are observed. The 
spatial density of young and intermediate age clusters de- 
creases toward the bulge, wh ich is consistent wi th a gas 
distribution in the M31 disk (|Nieten et al.ll2006D . 

To explain the M31 spiral structure and features of the 
prominent star-forming ring at i?M3i 10 kpc, two sce- 
narios of the M32 galaxy having passed through the M31 
disk were recently prop osed: (1) passage tro ugh the disk 
occurred ~20Myr ago () Gordon et al. I [20061 ) and created 
a "split" of the ring around position [6', 16'] (see Fig- 
ure 7; NGC 206 is located at the ring's "branching" point 
[15', 17']); (2) pass age trou gh the M31 center occurred 
-210 Myr ago (Bl ock et al.l l2006). The star cluster pop- 
ulation data could provide a reference to disentangle the 
more preferable scenario. 

True locations of objects in the M31 disk are subject to 
a strong projection effect due to an unfavorable inclina- 
tion angle of —75°. The projection of the galaxy and its 
rotation are such that the southwestern pa rt of the disk 
approaches the observer (|Hendersoiilll979D: the corota- 
tion radius l ies beyond i?M3i ~ 20 kpc (jEfremovl [19801 : 
IIvanov|[l985f ): and the rotation period at i?M3i ~ 10 kpc 
distance is of ~250 Myr. Assuming that all clusters in our 
sample are residing in the M31 disk, we made a deprojec- 
tion of cluster locations to face-on view, and referring to 
the disk rotation curve (jCarignan et al.|[2006l ). computed 
object positions as a function of look-back time. 

Judging by eye, the clumps of young clusters, which are 
located at [9', 5'] and [17', 10'], formed two compact spa- 
tial configurations — 80± 20 Myr ago. This age coincides 
well with the cluster age peak (see Section |4?2|) and sug- 
gests a possible evolutionary connection. There was no 
other compact spatial cluster configuration noted in the 
look-back time range from ~20Myr to ~150Myr, which 
is valid for this kind of study based on our cluster sam- 
ple. The range of the age distribution peak, estimated 
from the SSP model fitting, is of -50 - 120 Myr. How- 
ever, there is an age-metallicity degeneracy observed at 
the age of —100 Myr, which could shghtly bias the age 
determination. Therefore, we note this interesting coin- 
cidence of cluster age peak and compact spatial cluster 
configurations. To discuss this finding more reliably, a 
global kinematic study of the M31 disk cluster popula- 
tion and a detailed dynamic model of the M31 and M32 
encou nter are needed. Note however, that lCaldwell et aT] 
(|2009f ) did not find evidence for a peak in star formation 
between 20 Myr and 1 Gyr. 

4.7. Cluster Formation and Disruption 

If star clusters are being formed continuously at a con- 
stant rate with a power-law cluster initial mass func- 
tion, their number increases with time in the mass ver- 
sus age diagram as described by iBoutloukos fc Lamerd 
(2003 0. Consequently, the mass of the most massive 
cluster in the older log age bins rises due to a statisti- 
cal effect of sample size. However, in our sample, we 
see a lack of massive log(m/m0) > 4 clusters older 
than —100 Myr (Figure 4), probably, indicating a non- 
constant cluster formation history. Therefore, using the 
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observed mass and age distributions of star clusters we 
attempt to derive their formation and disruption rates in 
the studied field by adopting the technique developed by 
iBoutloukos fc LamersI (|2003l ). 

For a simple estimate, first we selected clusters more 
massive than log(TO/mQ) ^ 3.5. There are A^i '~ 30 
and N2 ^ 70 objects in the age ranges Aii = 30 < 
t < 100 Myr and Ais = 100 Myr < t < 1 Gyr, re- 
spectively. If a cluster formation rate (CFR) is con- 
stant and there is no cluster disruption, then a ratio of 
N2/N1 — At2/Ati « 13 is expected. The actual ratio is 
N2/N1 2, indicating that CFR is not constant and/or 
significant cluster disruption took place during the last 
Gyr in the M31 disk. 

The number of clusters per linear age interval, diV/dt, 
plotted versus age in Figure 5(a) can be described by two 
power-law slopes, intersecting at t ~ 300 Myr. The slope 
values were taken from the f ading/disruption models by 
IBoutloukos fc LamersI (|2003l Equation (14): a = 2.0, 
7 = 0.62, ( = 0.69) and only vertical adjustment to the 
observed distribution was applied. In the age domain of 
t < 300 Myr the slope line, expected for a cluster popula- 
tion forming continuously and fading in the course of evo- 
lution below the detection limit, provides a rather good 
match to the observed data. Around t ~ 70 Myr, there 
is a sudden increase in the number of clusters by a factor 
of ~2, which could be attributed to the increased CFR 
at that epoch. In the age domain older than ^300 Myr, 
the observed distribution can be well described by the 
slope expected for a cluster disruption at some arbitrary 
time. An increase in the cluster number at t 10 Gyr is 
attributed to the GC candidates, prominent in our sam- 
ple. 

The number of clusters per linear mass interval, 
diV/dm, is displayed versus mass in Figure 5(b). In the 
mass range of Iog(r7i/m0) < 3.3 the incompleteness of 
our cluster sample is obvious. Due to incompleteness, the 
mass distribution of our cluster sample cannot be used 
to constrain the typical lifetime of star clusters since the 
fading line cannot be adjusted properly. Note that the 
distribution of clusters in between ~ 100 Myr and ^--^3 Gyr 
and in the mass range of log(m/TOQ) > 4.5 is steeper 
than that of the whole sample, which includes the old 
GC candidates. 

Using the IBoutloukos fc LamersI (|2003D technique and 
relying on the power-law slope line intersection point in 
Figure 5(a), we estimate a conservative value of typi- 
cal lifetime of log(m/m0) = 4 mass cluster to be of 
^4^'^ ~ 300 Myr in the M31 survey field. For compari- 
so n, the following valu es in other galaxies were derived 
bv ILamers et al.l ()2005D : M51 central region— -70 Myr; 

M33 630 Myr; MW solar neighborhood 560 Myr; 

and SMC 8 Gyr. 

The cluster disruption analysis gives the disruption 
time of '^500 Myr if only cluster age distribution of our 
sample is used. If the mass distribution is taken into ac- 
count (although it is difficult to define intersection point 
of fading and disruption lines in Figure 5(b)), the clus- 
ter disruption time can be estimated to be as small as 
~ 100 Myr since it is sensitive to the "intersection" mass. 
Therefore, we choose tf'^ ~ 300 Myr as a conservative 
value. We note that a similar inconsistency between 
age and mass distribution analysis was noticed recently 



for LMC star clusters bv iParmentier fc de GriisI (|2008f) . 
However, they found an opposite effect for the LMC 
sample — the cluster mass distribution analysis tends to 
yield higher cluster disruption timescale than the age dis- 
tribution analysis. Such discrepancy might be related 
to the assumptions applied in the Boutloukos fc Lameri 
(2003) analysis used in our study: (1) a constant cluster 
formation rate; (2) a power-law cluster initial mass func- 
tion, which should be of the Schechter's function type as 
s hown in S ection 14.31 

ILamers et al. (2005) found anticorrelation between the 
hfetime, tf^^, of log(m/m0) = 4 mass cluster and the 
ambient density of its environme nt, Pambi in th e galax- 
ies. Following the reasoning by ILamers et al.l ([2005), 
we attempt to estimate Pamb in the M31 disk. We 
took the structural par ameters of the M31 disk from 
the "best- fit model" of IGeehan et all (|2006D : a disk 
scale-length, Rd = 5.4 kpc; a central surface density. 
So = 4.6 X 10^ 7710 kpc^^. Therefore, the average disk 
surface density at i?M3i ^ 10 kpc, i.e., at the center of 
our survey field S |=So exp(— i?M3i/-Rd) ~ 70 7770 pc~^. 

[van der KruitI (j2002l ) has shown that for spiral galaxies 
the vertical scal e heigh t is of hz ^ 0.15 i?^; this implies 
/jstars _ 800 pc. IBraunI ([l99l found /7f" - 350 pc for a 
gas distribution at a galactocentric distance of Rm31 = 
10 kpc. The scale height for dust di stribution is assumed 
to b e of fe^""* ~ 100 a nd ~1 50pc bv lHatano et"al] (llQQTh 
and ISemionov et alj (|2003D . respectively. Depending on 
the adopted hz, the estimated ambient density in our 
survey field, Pamb ~ 0.5Tth~^, is in the range of Pamb ~ 
0.05 ... 0.35 70,0 pc^'^. This density is much lower than 
the typical half-mass density ph ~ 10^ 7)7,0 pc~'^ of star 
clusters, shown in Figure 3(b)^ 

Referring to Figure 4 in ILamers et al.l (j2005f ). we 
see that values of tf^^ ~ 300 Myr and Pamb 
0.05 . . . 0.35 7770 pc~'^ derived here, place the studied field 
of the M31 disk in the position between MW and M33 
and M51, however, slight ly below the predic ted theoret- 
ical tf^ = if "(pamb) line. (Cieles et all (|2006[ ) have mod- 
eled encounters of star clusters with molecular clouds 
to explain short disruption time for star clusters in the 
M51 central region. Our survey field in M31 is cen- 
tere d on the star-form ing ring possessing high-density 
gas (jNieten et al.ll20d6l ). therefore, enhanced star cluster 
disruption in respect to MW and M33 co uld apparently 
take place. Although iRosolowskvl (|2007l and references 
therein) discuss similarity of the molecular cloud prop- 
erties in M31 and MW, enhanced star cluster disruption 
should apparently take place in the M31 molecular gas 
"ring" in respect to that of the solar neighborhood. The 
whole field survey of the M31 disk clusters and more de- 
tailed models of cluster and cloud kinematics, which are 
needed to estimate encounter cross sections in the disk, 
should be considered in future. 

Finally, we attempt to estimate the CFR in the M31 
disk. We select = 51 clusters with age At = 30 Myr < 
t < 130 Myr and mass 3.3 < log(777/7770) < 4.5 (see Fig- 
ure 4) . This subsample should have accurately estimated 
parameters and be relatively free of selection and cluster 
disruption effects. The total stellar mass of these clusters 
is Am ~ 3 X 10^ 77i0. Assuming t he cons tant power-law 
cluster mass function (index -2.0: lGielejl2009 ) down to 
log(777/77i0) = 2, mass correction factor of ~2.5 was de- 
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duced. Therefore, the approximate rate of star formation 
in the clusters was CFRgeid = Am/ At ^ 0.008 WQyr"^. 
The deprojected area of survey field (Figure 7) is of 
~70kpc^, making ~15% of the galaxy disk. From this, 
we estimate the lower limit of the average CFR at the 
epoch of ~100 Myr to be of CFRmsi ~ 0.05 m© yr~^ over 
the M31 disk. This can be compared t o the present- 
day star formation rate of 0.4moyr~^ (jBarmbv et al.l 
|2006[ ). indicating that >10% of formed stars could re- 
main "locked" in star clusters during ^100 Myr. 

For comparison, Williams (,2003) found the mean star 
formation rate for the disk ~lm0yr~^ over the last 
60 Myr. Based on stellar population analysis, he suggests 
that the lowest star formation rate occurred '^25 Myr 
ago being ~2 times lower compared with the epoch from 
~50Myr to ~250Myr. Note, a good coincidence of star 
formation history in M31 deduced by iWilliami ()2003f ) 
and one inferred from our cluster population data (Fig- 
ure 5(a)), which suggests an episode of a double increase 
in cluster formation ^^70 Myr ago. Also, iKodaira et al.l 
(jT999) have found that the star formation rate in the 
field around OB association A24 has decreased by ~2.5 
times from the epoch at '-~^1 Gyr to the present day. 

Although clusters in the intermediate mass range 
(3.5 < log(m/mo) < 4.5) are not well known in the MW, 
the compact clusters in M31 might represent a new class 
of star clusters of intermediate age 30 Myr ^ t < 3 Gyr, 
linking the age of GC formation and the present. It could 
be suspected that this class of disk-population clusters 
is not observed in the MW because of high extinction 
throug h the Galac t ic pla ne. This might be supported 
by the iClark et al.l (|2005 l) discovery that a Galactic star 
cluster Westerlund 1 is a Super Star Cluster with a mass 
of up to log(m/m0) ~ 5, age of ~3Myr, and radius 
of 0.6 pc, which suffers strong interstellar extinction of 
Ay ^ 11.5 mag being at the distance of ~5kpc. There- 
fore, a more numerous massive compact star cluster pop- 
ulation also could be hidden by dust clouds in the M31 
disk. 

5. SUMMARY 

The sample of compact star clusters in the present ho- 
mogeneous photometric survey of the 17.5' x 28.5' field in 
the southwestern part of the M31 disk is apparently over- 
lapping with the clusters detected in the patchy HST 
fields, in terms of their mass 3.0 < log(m/mo) < 4.5, 
which is in-between of the mass of typical DCs and clas- 
sical GCs, well known in the MW. The structural prop- 
erties of the compact clusters are similar to those of GCs, 
except for minor irregularities, and they resemble those 
of massive clusters in the LMC/SMC, including blue clus- 
ters of globular appearance. The metallicity of sample 
clusters is rather constant over a wide range of galacto- 
centric distances, with a typical value scattering around 
[M/H] 0.4 dex. 

The global follow-up survey of the whole M31 disk clus- 
ter population, based on photometry and size measure- 
ment, is of great interest and importance to comprehend 
the evolution of galaxy disks. Supplemented with spec- 
troscopic and higher image resolution studies, as well as 
modeling of cluster stochastic effects, it would provide 
a guiding panoramic view and better constraints for the 
M31 and M32 interaction models. 
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APPENDIX 
OBJECTS FOR DETAILED STUDY 

Here, we present several remarkable objects, which 
could serve as targets for future spectroscopy and high- 
resolution imaging. We briefly discuss their properties 
referring to the derived parameters, multiband and HST 
images. We provide cross-referencing with The Revised 
Bologna Catalog of M31 Globu lar Clusters and candi- 
dates compiled bv ,Galleti et al.l (12007) in brackets when 
available. 

KW044 (B325) is located at the position of [16',2'] in 
Figure 7 and is resolved in Suprime-Cam images. The 
structural model fit resulted in a large half-light radius, 
rfi ^ 11 pc, and an extremely small ratio of King model 
parameters, rt/rc ^ 3. Its age of '^60Myr, mass of 
log(TO/mQ) ~ 4.6, and radial velocity of —560km s~^, 
compatible with gas disk velocity at this position, im- 
ply that KW044 is a representative of the young disk 
cluster population. However, another solution for evolu- 
tionary parameters (a low extinction case) also exists — 
the cluster could be an old object of low concentra- 
tion and low metallicity, res embling the extended ones 
in M 31 (|Mackev et al.l I2OO6D and M33 (jStonkute et all 
|2008() , and suggesting its loca tion in the M31 halo. The 
age of ^630 Myr is derived bv lCaldwell et all (|2009( ). 

KW072 (VdBO) is located at the position of [15', 7'] in 
Figure 7. A visual impression that it is surrounded by an 
enhanced surface number density of stars (ejected from 
cluster) is confirmed by the structural model fit, indicat- 
ing its extended nature — structural model parameters: 
ratio rt/rc ^ 300 (King), n ^ 1 (EFF), and rh ^ 5.5 pc. 
Its young age of ^20 Myr is confirmed by a strong ul- 
traviolet [GALEX) flux. A mass, log(TO/m0) ~ 4.6, of 
the cluster, which formed recently, is remarkably high 
comparing with the MW t ypical OCs' m ass. This clus- 
ter was studied in detail bv lPerina et al.l ([2009) based on 
photometry of individual stars from HST images. The 
parameters derived are in a very good agreement with 
values derived in the present study based on its inte- 
grated properties. 

KW141 (BOllD) is located at the position of [11', 15'] 
in Figure 7. It represents a population of young, 
~10Myr, intermediate-mass, log(m/m0) ~ 3.8, clusters. 
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Fig. 8. — HST images. Panels show: (a) KW141 — an example 
of stochastic appearance of red bright stars (color coding of HST 
filters in the electronic edition: red — F8I4W, green — F555W, and 
blue—F439W}; (b) KW249— an object with the largest half-light 
radius, r^i ^ 14 pc, in our sample (F606W filter, overplotted with 
ellipse indicating r^, derived via elliptical King model fit). Note 
different image scales in panels (a) and (b). North is up, east is 
left. 

However, stochastically appearing four bright red stars 
enclose the blueish central part of this object suggesting 
an older age. Its HST WFPC2 color image is displayed 
in Figure 8(a). The age and mass of this "stochast ic clus- 
ter ex ample" are in good agreement with Caldwe ll et al.l 
(|2009f) : ~30Myr and log(TO/m0) ~ 3.8. 

KW249 — the most controversial object, located close 
to the M31 bulge at the position of [3', 29'] in Figure 7. 
Its young age of ~10Myr and a mass of log(TO/m0) ^ 4.3 
make it incompatible with the structural parameter best- 
fit model indicating very high ratio of rt jvc > 500 (King 
model) and large half-light radius of ~ 14 pc (the most 
extended object in our cluster sample). It is associated 
with 24 /im (Spitzer) and HI emission, probably, still be- 
ing an embedded cluster with E{B — V) ^ 1.5, which 
could explain the absence of Ha and ultraviolet emis- 
sion. It was imaged with the HST WFPC2 (F606W 
filter), displayed in Figure 8(b) (fitted elliptical models 
provide the ellipticity of ~0.75 — a ratio of minor to ma- 
jor axis). We suggest, that it could be either: (1) an 
intermediate-mass cluster, which expanded after a rapid 
gas removal phase (Baumgardt & Kroupa 2007), or (2) 
a background galaxy, having emission in Kg (2MASS) 
attributed to older stellar populations. 



KW271 & KW273 (G099*)—& double cluster can- 
didate, located at the position of [2', 18'] in Figure 7. 
KW271 is centered on a strong 24 /xm emission source, 
its parameters are: t ^ lOMyr, log(m/m0) ~ 3.7, 
E{B - F) ~ 0.9, and rh - 3.6 pc. Contrarily, KW273 
has no 24 /im emission, is much older {t ~ 800 Myr) , 
more massive (log(m/m0) ~ 4.5), has a low color ex- 
cess {E{B — V)^ 0.1), and large half-hght radius (r/j ~ 
4.3 pc). The structural parameters of two "adhered" 
cluster candidates were derived by masking KW271 and 
fitting the brighter KW273. After subtracting the best-fit 
model of KW273 from the original image, KW271 was 
fitted. In total, seven iterations of the object subtrac- 
tion were performed, leading to fine residual images. We 
suggest two possible explanations of the nature of these 
objects: (1) they make a double cluster, similar to those 
detected in LMC/SMC ijCarvalho et al.1 12008'). with one 
component just emerging from the embedded phase. 
Close image inspection reveals an extended "halo" en- 
shrouding both these objects, however, the derived age 
difference does not support this case; (2) both clus- 
ters (they could be a cluster and a galaxy — 24 /j,m emis- 
sion source) are subject to projection, however, proba- 
bility of such close coincidence is low. We note, that 
ICaldwell et all (.2009 ) classify KW271 as a background 
galaxy based on the measured redshift z ~ 0.15, and 
for KW273 they provide age of ~400Myr and mass 
log(m/m0) 4.0 in reasonable agreement with our val- 
ues. 

KW279 (B205D) is located at a large projected galac- 
tocentric distance, thus far from the star-forming ring, 
at the position of [2', 5'] in Figure 7. It is associated with 
a strong 24 /^m emission, which supports its young age 
of '^lOMyr. A low-mass, log(m/m0) ~ 3.2, substan- 
tial color excess, E{B — V)^ 0.6, and a moderate size, 
rh 2pc, were deduced assuming it being a star clus- 
ter. However, its asymmetric shape with a blueish "tail," 
pointing at a neighboring 24 /im emission source, implies 
that it cou ld be an in t eracti ng starburst galaxy. We note, 
that Cald well et all (|2009D classify KW279 as a back- 
ground galaxy based on the measured redshift z ^ 0.1. 
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TABLE 1 

Evolutionary and Structural Parameters of High- Probability Star 
Cluster Candidates in the Southwest Part of the M31 Disk 



ID 


aj20oo 


5j2000 




log(t/Myr)b 


log(m,/mQ 


[M/H]<1 


E{B - V)" 




KW003 


10.04519 


40.89759 


-6.0 


1.89 


3.42 


-0.5 


0.12 


2.0 


KW004 


10.04568 


40.60326 


-6.6 


1.67 


3.48 


-0.9 


0.10 


0.7 


KW006 


10.05430 


40.60535 


-6.7 


1.79 


3.66 


-0.2 


0.21 


0.7 


KW007 


10.05468 


40.69164 


-6.1 


2.04 


3.56 


-0.4 


0.52 


0.9 


KW008 


10.05671 


40.63810 


-6.1 


1.52 


3.16 


-0.8 


0.30 


0.6 


KW009 


10.05700 


40.76773 


-5.2 


1.36 


2.79 


-0.1 


0.34 


0.8 


KWOlO 


10.05768 


40.67549 


-7.3 


2.07 


4.01 


-0.5 


0.70 


5.1 


KWOll 


10.05924 


40.65584 


-6.6 


1.88 


3.66 


-0.8 


0.16 


1.4 


KW013 


10.06067 


40.62244 


-6.6 


1.78 


3.61 


-0.4 


0.02 


1.5 


KW014 


10.06177 


40.77234 


-6.3 


1.42 


3.26 


-0.1 


0.50 


2.2 


KW015 


10.06203 


40.75771 


-4.9 


2.85 


3.70 


0.1 


0.05 


1.9 


KW017 


10.06453 


40.66658 


-6.9 


1.72 


3.64 


-0.7 


0.22 


1.0 


KW019 


10.06504 


40.58820 


-5.2 


2.56 


3.54 


-0.3 


0.13 


1.3 


KW020 


10.06519 


40.59857 


-5.1 


3.39 


4.07 


-1.1 


0.09 


2.6 


KW022 


10.07194 


40.65138 


-7.0 


1.74 


3.72 


-0.4 


0.03 


3.0 


KW023 


10.07237 


40.79158 


-6.8 


1.60 


3.55 


-0.5 


0.41 


2.7 


KW024 


10.07249 


40.54608 


-5.1 


1.59 


2.92 


-0.1 


0.20 


1.2 


KW025 


10.07289 


40.58077 


-7.0 


3.71 


5.02 


-1.9 


0.52 


8.7 


KW026 


10.07313 


40.65559 


-6.6 


1.83 


3.63 


-0.8 


0.05 


1.5 


KW027 


10.07362 


40.55231 


-6.1 


1.89 


3.51 


0.0 


0.25 


2.6 


KW028 


10.07611 


40.54573 


-7.0 
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^J2000 
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KW039 


10.08462 


40.73291 


-9.3 


3.82 


6.15 


-0.6 


0.15 


1.2 


KW040 


10.08659 


40.55609 


-6.5 


1.97 


3.66 


-0.5 


0.15 


1.0 


KW041 


10.08782 


40.62758 


-4.7 


1.96 


2.92 


-0.4 


0.03 


1.3 


KW042 


10.08929 


40.62074 


-5.9 


2.69 


3.92 


-0.6 


0.21 


2.2 


KW043 


10.09325 


40.63796 


-6.3 


1.40 


3.24 


-0.2 


0.35 


0.7 


KW044 


10.09623 


40.51320 


-9.2 


1.77 


4.59 


-0.8 


0.60 


10.0 


KW045 


10.09674 


40.62105 


-5.9 


2.27 


3.68 


-0.1 


0.22 


1.1 


KW046 


10.09793 


40.64923 


-6.3 


1.45 


3.19 


-0.8 


0.18 


0.7 


KW048 


10.10042 


40.60629 


-8.4 


1.30 


4.02 


-0.2 


0.19 


0.8 


KW049 


10.10088 


40.53335 


-6.0 


1.92 


3.44 


-0.4 


0.03 


1.5 


KW050 


10.10110 


40.68546 


-5.4 


1.72 


3.07 


-0.4 


0.03 


1.7 


KW051 


10.10121 


40.61724 


-4.5 


2.26 


3.13 


0.1 


0.02 


0.9 


KW053 


10.10341 


40.81300 


-7.6 


2.40 


4.40 


-0.5 


0.37 


7.5 


KW054 


10.10363 


40.81692 


-7.8 


2.01 


4.20 


-0.4 


0.69 


4.6 


KW055 


10.10418 


40.55488 


-5.5 


1.91 


3.24 


-0.4 


0.02 


2.4 


KW056 


10.10534 


40.81352 


-6.9 


2.41 


4.15 


-0.4 


0.38 


2.8 


KW057 


10.10801 


40.62815 


-9.8 


0.84 


3.97 


-0.5 


1.15 


6.8 


KW058 


10.10906 


40.75886 


-6.1 


1.69 


3.34 


-0.2 


0.47 


1.0 


KW059 


10.11086 


40.79259 


-7.4 


2.11 


4.05 


-0.8 


0.46 


1.8 


KW061 


10.11223 


40.53262 


-7.3 


1.99 


3.99 


-0.4 


0.63 


2.6 


KW062 


10.11367 


40.75677 


-6.7 


1.98 


3.73 


-0.6 


0.12 


1.4 


KW064 


10.11449 


40.66775 


-5.9 


1.74 


3.28 


-0.5 


0.24 


1.1 


KW065 


10.11469 


40.51168 


-4.9 


2.46 


3.37 


0.0 


0.03 


3.2 
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10.11525 


40.52838 


-5.6 


1.48 


2.99 
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0.05 
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10.11554 


40.77210 


-7.7 


1.04 


3.48 
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40.54411 


-5.0 


2.71 


3.61 
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40.77968 


-4.7 


1.49 


2.61 
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0.04 
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KW075 


10.12328 


40.54568 


-5.2 


3.12 


3.84 


-1.2 


0.06 


1.5 


KW076 


10.12387 


40.63842 


-5.8 


2.58 


3.81 


-0.4 


0.30 


1.7 


KW077 


10.12448 


40.51532 


-5.7 


2.55 


3.74 


-0.2 


0.27 


3.2 


KW078 


10.12708 


40.75818 


-8.5 


1.18 


3.97 


-0.5 


0.03 


1.0 


KW079 


10.12764 


40.74831 


-7.1 


1.11 


3.54 


0.5 


0.03 


1.7 


KW081 


10.12861 


40.83648 


-6.7 


2.40 


4.04 


-0.5 


0.50 


4.3 


KW082 


10.12872 


40.78416 


-5.6 


2.03 


3.38 


-0.3 
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10.13069 


40.75239 


-5.1 


0.95 


2.52 
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10.13216 


40.51760 


-4.8 


2.30 


3.21 
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KW086 


10.13340 


40.74484 


-5.9 


0.70 


2.44 


-0.0 


0.09 
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KW087 


10.13404 


40.77417 


-6.0 


3.76 


4.64 


-1.8 


0.26 


7.4 


KW089 


10.13572 


40.83714 


-6.3 


2.29 


3.83 


-0.3 


0.07 


2.2 


KW094 


10.14379 


40.80833 


-6.7 


1.11 


3.19 


0.1 


0.65 


0.7 


KW096 


10.14847 


40.93010 


-6.8 


1.18 


3.30 


-0.5 


0.83 


1.0 


KW097 


10.14988 


40.65192 


-6.9 


1.81 


3.73 


-0.9 


0.47 


1.1 


KW099 


10.15130 


40.84943 


-6.1 


3.44 


4.52 


-1.2 


0.18 


3.1 


KWlOO 


10.15220 


40.67090 


-8.0 


1.94 


4.27 


-0.4 


0.51 


5.6 


KWlOl 


10.15475 


40.55608 


-6.2 


2.82 


4.08 


-1.0 


0.42 


1.5 


KW102 


10.15519 


40.65407 


-7.5 


1.76 


3.92 


-1.0 


0.55 


1.3 


KW103 


10.15562 


40.81270 


-6.4 


2.23 


3.84 


-0.3 


0.10 


3.7 
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10.15837 


40.77468 


-4.8 


2.00 


3.03 


-0.3 


0.06 


1.5 
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10.16050 


40.74816 


-5.5 


3.11 


3.96 


-0.9 


0.26 


2.7 


KW107 


10.16147 


40.52894 


-4.5 


2.67 


3.36 


-0.1 


0.08 


1.1 


KWllO 


10.16421 


40.80204 


-4.9 


3.02 


3.74 


-0.4 


0.05 


3.3 


KW112 


10.17363 


40.64115 


-8.2 


3.88 


5.71 


-1.2 


0.51 


6.0 


KW113 


10.17612 


40.60131 


-6.5 


3.38 


4.67 


-0.7 


0.06 


6.5 


KW114 


10.17671 


40.76094 


-8.7 


2.00 


4.56 


-0.3 


1.16 


2.9 


KW115 


10.17808 


40.83946 


-5.3 


2.05 


3.27 


-0.3 


0.06 


2.7 


KW119 


10.18446 


40.56793 


-7.2 


1.38 


3.55 


-0.6 


0.48 


2.8 


KW120 


10.18520 


40.62048 


-8.8 


1.94 


4.57 


-0.7 


1.07 


2.6 


KW121 


10.18651 


40.50836 


-7.2 


1.99 


3.93 


-0.6 


0.65 


3.4 


KW122 


10.18702 


40.88564 


-6.5 


2.58 


4.06 


-0.6 


0.41 


1.6 


KW123 


10.19208 


40.95922 


-5.4 


2.30 


3.51 


0.1 


0.03 


2.3 


KW124 


10.19336 


40.86136 


-6.7 


1.81 


3.65 


-0.8 


0.30 


2.5 


KW125 


10.19582 


40.68280 


-5.7 


3.42 


4.38 


-0.6 


0.04 


3.2 


KW126 


10.19808 


40.66878 


-8.9 


0.70 


3.64 


0.1 


1.37 


3.2 


KW127 


10.19883 


40.63281 


-6.2 


3.16 


4.31 


-0.4 


0.24 


5.9 


KW129 


10.19919 


40.92626 


-5.8 


2.38 


3.71 


-0.0 


0.05 


4.5 


KW130 


10.20018 


40.51726 


-5.7 


3.32 


4.27 


-0.8 


0.07 


2.9 


KW131 


10.20031 


40.67784 


-6.9 


2.31 


4.04 


-0.6 


0.77 


2.6 


KW132 


10.20148 


40.86623 


-7.0 


1.34 


3.51 


-0.1 


0.23 


1.3 


KW133 


10.20148 


40.58504 


-6.8 


2.31 


3.99 


-0.5 


0.10 


1.9 


KW134 


10.20233 


40.96013 


-6.7 


2.50 


4.12 


-0.2 


0.33 


2.1 


KW135 


10.20362 


40.56563 


-6.5 


1.83 


3.64 


0.0 


0.19 


6.5 


KW136 


10.20492 


40.64867 


-5.5 


2.93 


3.90 


-1.0 


0.07 


1.1 


KW137 


10.20583 


40.69227 


-7.5 


3.36 


5.01 


-1.0 


0.18 


4.2 
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FM /Hid 
[ivi/nj 


E{^B — y)^ 
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KW139 


10.21178 


40.67737 


-6.0 


3.20 


4.36 


-0.5 


0.39 


3.3 


KW140 


10.21457 


40.55768 


-7.7 


2.31 


4.33 


-0.6 


0.47 


6.7 


KW141 


10.21504 


40.73508 


-8.7 


0.84 


3.81 


0.6 


0.65 


0.8 


KW142 


10.21700 


40.58377 


-5.8 


1.85 


3.30 


-0.5 


0.32 


0.7 


KW143 


10.21761 


40.53484 


-4.9 


2.66 


3.49 


-0.5 


0.06 


3.1 


KW144 


10.21773 


40.97823 


-6.6 


2.78 


4.21 


-1.0 


0.33 


1.8 


KW145 


10.21777 


40.89900 


-6.5 


3.62 


5.04 


-0.1 


0.07 


8.8 


KW146 


10.22036 


40.86882 


-6.2 


3.20 


4.33 


-1.3 


0.35 


4.2 


KW147 


10.22061 


40.58877 


-7.1 


3.75 


5.11 


-1.9 


0.34 


5.8 


KW149 


10.22197 


40.73255 


-6.5 


0.70 


2.66 


-0.0 


0.62 


1.0 


KW150 


10.22235 


40.82742 


-6.1 


3.81 


4.87 


-0.6 


0.42 


3.8 


KW151 


10.22301 


40.71887 


-5.9 


2.37 


3.70 


-0.3 


0.51 


1.2 


KW152 


10.22375 


40.61417 


-5.4 


3.97 


4.62 


-1.4 


0.24 


1.1 


KW153 


10.22448 


40.671,37 


-4.6 


2.75 


3.49 


-0.1 


0.06 


1.2 


KW154 


10.22840 


40.59210 


-5.5 


2.64 


3.68 


-0.8 


0.14 


1.6 


KW155 


10.22846 


40.73891 


-5.4 


3.37 


4.19 


-0.8 


0.09 


2.0 


KW156 


10.23071 


40.62244 


-5.9 


1.90 


3.47 


0.2 


0.49 


2.1 


KW158 


10.23358 


40.70630 


-7.0 


2.20 


3.99 


-0.6 


0.62 


1.6 


KW159 


10.23475 


40.57826 


-5.3 


2.34 


3.42 


-0.6 


0.21 


1.2 


KW160 


10.23602 


40.57370 


-5.9 


2.48 


3.77 


-0.6 


0.06 


3.6 


KW161 


10.23726 


40.60918 


-7.3 


2.45 


4.29 


-0.5 


0.83 


2.2 


KW164 


10.23953 


40.74093 


-7.8 


4.01 


5.58 


-1.9 


0.59 


1.8 


KW165 


10.23976 


40.54874 


-6.2 


1.58 


3.32 


-0.2 


0.25 


1.9 


KW166 


10.24050 


40.57272 


-5.3 


1.43 


2.84 


-1.0 


0.06 


0.7 


KW168 


10.24321 


40.89079 


-9.3 


0.70 


3.81 


0.4 


1.49 


2.3 


KW169 


10.24352 


40.60981 


-6.9 


1.50 


3.58 


-0.0 


0.39 


0.8 


KW170 


10.24504 


40.57345 


-6.2 


1.85 


3.59 


0.3 


0.27 


0.9 


KW171 


10.24531 


40.59666 


-10.6 


4.00 


6.72 


-1.6 


0.14 


3.7 


KW172 


10.24540 


40.71805 


-6.4 


3.83 


5.23 


-0.0 


0.51 


4.9 


KW173 


10.24663 


40.58450 


-5.9 


1.86 


3.34 


-0.7 


0.26 


2.4 


KW174 


10.24728 


40.56506 


-5.3 


2.38 


3.48 


-0.4 


0.18 


1.1 


KW176 


10.25111 


40.57447 


-7.4 


0.90 


3.09 


-0.3 


0.81 


1.4 


KW177 


10.26062 


40.84083 


-5.0 


2.12 


3.14 


-0.5 


0.19 


1.0 


KW178 


10.26071 


40.80432 


-6.8 


4.02 


5.51 


-0.0 


0.34 


3.2 


KW180 


10.26142 


40.56423 


-6.1 


2.03 


3.56 


-0.5 


0.14 


1.8 


KW181 


10.26186 


40.69644 


-5.9 


2.84 


4.00 


-0.9 


0.36 


1.9 


KW182 


10.26203 


40.58293 


-6.5 


1.93 


3.65 


-0.4 


0.14 


1.2 


KW183 


10.26229 


40.94328 


-5.6 


4.09 


5.11 


-0.0 


0.08 


9.7 


KW184 


10.26513 


40.92902 


-6.3 


3.30 


4.45 


-1.4 


0.36 


5.8 


KW185 


10.26663 


40.53928 


-6.8 


2.26 


3.96 


-0.5 


0.43 


9.2 


KW186 


10.26836 


40.66409 


-5.2 


2.98 


3.84 


-0.5 


0.12 


1.0 


KW187 


10.26897 


40.57899 


-6.5 


1.23 


3.22 


-0.1 


0.20 


1.2 


KW188 


10.27108 


40.75779 


-6.1 


2.46 


3.85 


-0.4 


0.38 


2.1 


KW189 


10.27233 


40.85868 


-5.2 


2.23 


3.32 


-0.5 


0.07 


1.5 


KW190 


10.27249 


40.58317 


-6.8 


1.54 


3.60 


0.1 


0.21 


0.9 


KW191 


10.27438 


40.89986 


-6.0 


2.65 


3.90 


-0.8 


0.24 


4.3 


KW192 


10.27591 


40.61567 


-6.1 


1.82 


3.41 


-0.8 


0.35 


3.7 


KW193 


10.27601 


40.60383 


-5.3 


2.39 


3.51 


0.1 


0.04 


3.6 


KW194 


10.27857 


40.57474 


-6.9 


1.80 


3.70 


-0.8 


0.24 


2.8 


KW195 


10.28042 


40.90625 


-4.7 


3.02 


3.66 


-0.4 


0.04 


1.8 


KW196 


10.28306 


40.88364 


-6.9 


4.01 


5.32 


-0.8 


0.11 


2.7 


KW197 


10.28319 


40.81240 


-6.5 


1.89 


3.61 


-0.7 


0.45 


5.4 


KW198 


10.28443 


40.54785 


-7.5 


2.10 


4.16 


-0.4 


0.51 


4.5 


KW199 


10.28814 


40.59810 


-8.4 


3.84 


5.81 


-0.7 


0.09 


2.2 


KW200 


10.29161 


40.96981 


-7.0 


3.19 


4.66 


-0.9 


0.16 


5.7 


KW201 


10.29176 


40.53308 


-4.8 


2.50 


3.40 


-0.0 


0.04 


3.3 


KW203 


10.29775 


40.80557 


-6.7 


1.34 


3.44 


0.2 


0.29 


1.6 


KW204 


10.29872 


40.58854 


-5.3 


2.32 


3.42 


-0.4 


0.16 


1.4 


KW205 


10.29920 


40.73882 


-5.0 


3.02 


3.85 


0.1 


0.04 


4.5 


KW206 


10.29984 


40.96378 


-6.9 


2.22 


3.94 


-0.6 


0.62 


4.0 


KW207 


10.30204 


40.67383 


-6.9 


2.09 


3.88 


-0.6 


0.66 


2.9 


KW208 


10.30326 


40.57155 


-7.7 


1.75 


4.01 


-0.3 


0.24 


2.9 


KW209 


10.30461 


40.64758 


-4.9 


2.76 


3.60 


-0.0 


0.06 


2.6 


KW210 


10.30754 


40.56614 


-7.3 


2.24 


4.19 


-0.4 


0.06 


9.3 


KW211 


10.31052 


40.93089 


-9.3 


4.00 


6.19 


-1.6 


0.17 


2.2 


KW212 


10.31492 


40.81869 


-5.2 


0.84 


2.09 


-0.2 


0.09 


1.3 


KW213 


10.31781 


40.86897 


-5.0 


2.93 


3.69 


-0.9 


0.03 


2.1 


KW214 


10.31870 


40.98434 


-8.2 


4.10 


5.78 


-2.2 


0.19 


1.7 


KW215 


10.31888 


40.65087 


-6.5 


1.84 


3.57 


-0.8 


0.14 


1.4 


KW216 


10.31985 


40.80730 


-4.5 


2.53 


3.30 


-0.0 


0.05 


2.5 


KW217 


10.32373 


40.72907 


-5.8 


2.49 


3.75 


-0.5 


0.32 


1.4 


KW218 


10.32443 


40.72654 


-7.9 


0.70 


3.16 


-0.3 


1.09 


2.9 


KW219 


10.32559 


40.73372 


-6.7 


2.79 


4.28 


-0.5 


0.20 


4.0 


KW221 


10.32802 


40.95444 


-8.1 


3.93 


6.03 


0.1 


0.35 


1.7 


KW223 


10.32991 


40.75117 


-6.8 


3.98 


5.20 


-1.4 


0.26 


1.3 


KW224 


10.33554 


40.60722 


-4.3 


1.84 


2.70 


-0.4 


0.02 


1.1 
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TABLE 1 — Continued 



ID 


'^,12000 


s 


72000 








nVT /HI'* 
[ivi/ nj 


E[B — vy 


' h 


KW225 


10.33715 


40, 


.98462 


-7.2 


3.95 


5.60 


-0.1 


0.17 


1.7 


KW226 


10.33754 


40, 


.73732 


-6.0 


2.30 


3.71 


-0.4 


0.41 


1.4 


KW227 


10.33773 


40, 


,70366 


-6.4 


2.40 


3.92 


-0.4 


0.39 


2.0 


KW228 


10.34021 


40, 


.68455 


-6.2 


2.23 


3.72 


-0.5 


0.28 


1.6 


KW229 


10.34250 


40, 


.83295 


-6.5 


3.79 


4.89 


-1.9 


0.48 


4.8 


KW230 


10.34662 


40, 


,73745 


-4.4 


2.88 


3.47 


-0.4 


0.05 


1.4 


KW231 


10.34725 


40, 


.63251 


-5.4 


2.42 


3.54 


-0.1 


0.24 


3.4 


KW232 


10.35035 


40, 


,61306 


-6.8 


1.98 


3.79 


-0.5 


0.06 


8.2 


KW233 


10.35365 


40, 


.75562 


-5.0 


2.92 


3.85 


0.3 


0.07 


1.3 


KW234 


10.35418 


40, 


.87940 


-4.6 


3.04 


3.77 


0.3 


0.07 


2.2 


KW235 


10.35582 


40, 


,51488 


-6.5 


3.82 


5.07 


-0.5 


0.56 


2.3 


KW236 


10.35815 


40, 


.56083 


-5.5 


2.58 


3.69 


-0.4 


0.20 


3.8 


KW239 


10.36204 


40, 


,78528 


-6.6 


0.95 


3.03 


0.1 


0.61 


3.1 


KW240 


10.36244 


40, 


,69372 


-6.8 


2.06 


3.81 


-0.4 


0.31 


1.8 


KW242 


10.36512 


40, 


.80364 


-6.3 


2.02 


3.60 


-0.4 


0.48 


0.9 


KW243 


10.36531 


40, 


,67627 


-5.0 


3.01 


3.75 


-0.5 


0.04 


1.1 


KW244 


10.36715 


40, 


,89721 


-9.1 


3.86 


6.22 


-0.3 


0.04 


1.8 


KW245 


10.37280 


40, 


.75469 


-6.1 


2.59 


3.95 


-0.3 


0.14 


4.2 


KW246 


10.37317 


40, 


,83282 


-7.8 


1.88 


4.19 


0.1 


0.94 


3.2 


KW247 


10.37492 


40, 


,64060 


-5.7 


1.76 


3.25 


-0.2 


0.20 


2.8 


KW248 


10.37625 


40, 


.80908 


-5.1 


2.47 


3.49 


0.1 


0.04 


4.8 


KW249 


10.37817 


40, 


.96697 


-10.5 


0.70 


4.30 


0.4 


1.46 


13.5 


KW250 


10.37882 


40, 


,65230 


-5.7 


3.40 


4.26 


-1.6 


0.24 


2.9 


KW251 


10.37981 


40, 


.77018 


-5.7 


2.56 


3.79 


0.1 


0.03 


2.2 


KW252 


10.38217 


40, 


,67094 


-6.5 


1.49 


3.34 


-0.5 


0.69 


2.1 


KW253 


10.38282 


40, 


.74846 


-4.6 


2.46 


3.27 


0.0 


0.05 


2.8 


KW254 


10.38501 


40, 


.81791 


-5.0 


3.07 


3.86 


0.0 


0.06 


3.0 


KW255 


10.38605 


40, 


,79066 


-5.6 


2.41 


3.65 


-0.1 


0.32 


2.9 


KW256 


10.38768 


40, 


,65660 


-5.6 


1.89 


3.25 


-0.5 


0.28 


1.6 


KW257 


10.38851 


40, 


.74583 


-6.3 


2.27 


3.76 


-0.5 


0.55 


3.9 


KW258 


10.38910 


40, 


,92238 


-9.9 


0.70 


4.05 


0.4 


1.66 


3.4 


KW259 


10.38952 


40, 


,66234 


-4.7 


2.57 


3.39 


-0.1 


0.05 


1.8 


KW260 


10.38954 


40, 


.61910 


-5.4 


2.38 


3.57 


0.0 


0.03 


2.3 


KW261 


10.38992 


40, 


,82127 


-8.5 


0.70 


3.48 


-0.0 


1.39 


2.1 


KW262 


10.39174 


40, 


.63835 


-5.2 


2.73 


3.72 


-0.1 


0.05 


5.0 


KW263 


10.39175 


40, 


.70641 


-7.0 


2.39 


4.16 


-0.4 


0.68 


2.8 


KW264 


10.39327 


40, 


,87083 


-5.8 


2.39 


3.73 


-0.1 


0.26 


1.9 


KW265 


10.39493 


40, 


.87377 


-6.6 


1.46 


3.38 


-0.2 


0.59 


1.1 


KW266 


10.39585 


40, 


.76342 


-6.0 


3.17 


4.24 


-0.6 


0.06 


8.6 


KW267 


10.39828 


40, 


,92745 


-6.7 


4.03 


5.50 


0.0 


0.37 


2.9 


KW268 


10.39855 


40, 


.97107 


-5.4 


4.01 


4.62 


-1.9 


0.16 


3.5 


KW269 


10.39923 


40, 


.58348 


-7.4 


1.08 


3.43 


-0.9 


0.99 


1.7 


KW270 


10.40237 


40, 


,96336 


-5.4 


2.60 


3.70 


0.1 


0.05 


1.0 


KW271 


10.40308 


40, 


.78978 


-9.3 


0.84 


3.75 


-0.2 


0.92 


5.4 


KW272 


10.40315 


40, 


,84670 


-9.0 


1.73 


4.50 


-0.7 


1.14 


9.0 


KW273 


10.40363 


40, 


,79043 


-7.0 


2.93 


4.52 


-0.5 


0.06 


6.1 


KW274 


10.40468 


40, 


.94694 


-4.9 


3.97 


4.73 


0.1 


0.07 


3.5 


KW275 


10.40531 


40, 


,68044 


-5.0 


2.76 


3.67 


0.1 


0.04 


2.7 


KW276 


10.40555 


40, 


,88961 


-5.0 


2.41 


3.38 


-0.1 


0.06 


2.3 


KW277 


10.40689 


40, 


.81064 


-6.6 


2.30 


3.88 


-0.7 


0.71 


2.3 


KW278 


10.40818 


40, 


,79498 


-6.4 


2.48 


4.02 


-0.2 


0.32 


2.0 


KW279 


10.40857 


40, 


,56961 


-7.8 


0.90 


3.23 


-0.4 


0.55 


1.9 


KW280 


10.41025 


40, 


.97424 


-6.0 


3.01 


4.25 


0.1 


0.15 


3.8 


KW281 


10.41050 


40, 


,82679 


-5.7 


2.77 


3.98 


0.2 


0.04 


1.5 


KW283 


10.41167 


40, 


.79971 


-6.3 


3.81 


5.05 


-0.3 


0.37 


8.1 


KW284 


10.41180 


40, 


.68181 


-7.4 


2.02 


4.06 


-0.5 


0.14 


1.2 


KW285 


10.41559 


40, 


,94267 


-6.5 


3.92 


5.26 


-0.1 


0.31 


2.8 



Note. — ID iiumbGr, R.A. & Dccl. (J2000) coordinates of the photometric aperture center in the USNO-Bl.O catalog system (degrees) according 
to lNarbutis et all ) |200gD . 

^ Absolute magnitude, Mv , corrected for aperture correction. 

^ Age, t in Myr, log(i/Myr). 

^ Mass, m in solar mass units, log(m/mQ). 

^ Metallicity, [M/H]. 

^ Interstellar extinction in MW plus M31; color excess, E{B — V). 

^ Average half-light radius in pc, — (t^ + r^)/2, derived from the King and EFF models. 



